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ABSTRACT 


Multidisciplinary design, analysis, and optimization using a genetic algorithm is being 
developed at the National Aeronautics and Space Administration Dryden Flight Research 
Center (Edwards, California) to automate analysis and design process by leveraging 
existing tools to enable true multidisciplinary optimization in the preliminary design stage 
of subsonic, transonic, supersonic, and hypersonic aircraft. This is a promising technology, 
but faces many challenges in large-scale, real-world application. This report describes 
current approaches, recent results, and challenges for multidisciplinary design, analysis, 
and optimization as demonstrated by experience with the Ikhana fire pod design. 

NOMENCLATURE 


AIC 

aerodynamic influence coefficient 

DC 

DOT optimizer with continuous design variable 

DFRC 

Dryden Flight Research Center 

DOT 

design optimization tool 

GA 

genetic algorithm 

GC 

genetic optimizer with continuous design variable 

GD 

genetic optimizer with discrete design variable 

KEAS 

knots equivalent airspeed 

MDAO 

multidisciplinary design, analysis, and optimization 

NASA 

National Aeronautics and Space Administration 

Pi,P 2 

loading values 

V-g 

speed versus damping 

V-(jl) 

speed versus frequency 

X r x 2 , x 3 

design variables 


INTRODUCTION 

Supporting the Aeronautics Research Mission Directorate (ARMD) guidelines, the 
National Aeronautics and Space Administration (NASA) Dryden Flight Research Center 
(DFRC) (Edwards, California) is developing a multidisciplinary design, analysis, and 
optimization (MDAO) tool. This tool will leverage existing tools and practices, such as MSC 
Nastran (MSC Software Corporation, Santa Ana, California), ZAERO® (Zona Technology, 
Inc., Scottsdale, Arizona), and computational fluid dynamics (CFD) codes, and allow the 
easy integration and adoption of new state-of-the-art software. 

Optimization has made its way into many mainstream applications. For example, 
MSC Nastran has developed solution sequence 200 for design optimization (ref. 1 ), and 


MATLAB (The MathWorks, Natick, Massachusetts) has developed an optimization toolbox 
(ref. 2). Other applications, such as the ZAERO® aeroelastic panel code (ref. 3) and the 
CFL3D Navier-Stokes solver (ref. 4) do not include a built-in optimizer. 

Current commercial MDAO tools are codes that are limited to certain disciplines and 
are largely based on a single fidelity approach. Aircraft designs using nonlinear analyses 
such as transonic aeroelasticity are predominantly developed using the manual trial-and- 
error approach; this considerably slows the entire design process. The development of 
an MDAO tool to automate the existing manual trial-and-error approach will enable true 
multidisciplinary optimization early in the design process. 

The primary and long-term objective of the development of the MDAO tool is to 
generate a “central executive” capable of using disparate software packages in a cross- 
platform network environment so as to quickly perform optimization and design tasks in a 
cohesive streamlined manner. This object-oriented framework will integrate the analysis 
codes for multiple disciplines, instead of relying on one code to perform the analysis for 
all disciplines. Optimization can then take place within each individual tool, or in a loop 
between the executive and the tool, or both. Figure 1 shows a typical set of tools and their 
relation to the central executive. 

BACKGROUND AND ANALYSIS MODULES 

The heart of the MDAO framework is the central executive module. This is the module 
the user will use to choose input files, solution modules, and output files; to determine the 
status of executing jobs; and to select modules for the viewing and filtering of outputs. It is 
a graphical user interface (GUI) that provides a single point of control for applications that 
run on the user’s own computer, or for code that may reside on remote workstations or a 
computational cluster. At this stage, most of the code is written in FORTRAN on a UNIX 
platform. The GUI will be added later in the development process. 

An interesting aspect of the central executive application is user-selectable resolution. 
An engineer engaged in preliminary design will probably prefer quick analysis turnaround 
over absolute accuracy. In aeroelastic analysis, for example, this may mean using a 
simple strip theory. As the design effort proceeds, the designer may segue into using a 
more accurate aeroelastic panel code such as ZAERO®, and then perhaps an unsteady 
Euler- or Navier-Stokes-type of solver. The interface and process, however, to the extent 
practicable, should be the same for each solution type. The user also can easily plug 
in his own analyzer and script. This could be considered a form of high-level object- 
oriented solution. Only the interface is exposed to the user. The interface variables can 
be total weight, safety factors, natural frequencies, drag, noise level, flutter speed, flutter 
frequencies, et cetera. Internally, one solution module could be entirely different from 
another, but the solution module is hidden from the user. 
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Two external optimizers are included in the MDAO tool: design optimization tools 
(DOTs) (refs. 5, 6) based on a gradient-based algorithm; and the genetic algorithm (GA) 
(ref. 7). One reason to use an external optimizer is to use a different type of optimizer. 
NASTRAN, for example, uses a gradient-based approach to optimization (ref. 8); a 
drawback to this approach is the necessity to compute finite difference or analytical 
sensitivity values and perform the search, and there is often a need for a prior input with 
experience in defining the problem and search directions. 

The DOT is a commercial optimization code that can be used to solve a wide variety 
of nonlinear optimization problems. When the DOT requires the values of the objective 
and constraint functions corresponding to a proposed design, it returns control to the 
user’s program. The user’s program calls the DOT again to obtain the next design point; 
this process is repeated until the DOT returns a parameter to indicate that the optimum 
objective function is reached. 

Genetic algorithms do not require gradient calculations and can be started with 
random seeds, eliminating some of the need for user input and allowing for solutions that 
may not be readily apparent to even experienced designers (ref. 9). In the case of multiple 
local minima problems, genetic algorithms are able to find the global optimum results 
while gradient-based algorithms may converge to the local optimum value. 

Multidisciplinary Design, Analysis, and Optimization Framework 

Currently, the central executive MDAO framework can handle structural optimization 
problems. The framework of the MDAO is object oriented: users can either use the built- 
in pre- and post-processor to convert design variables to structural parameters and 
generate objective functions, or use their own analyzer for the optimization analysis. The 
MDAO framework process is presented in the flowchart in figure 2. The MDAO tool has 
been tested and validated in several test cases and applied to model update problems 
(ref. 10). Basically, the current MDAO framework developed at NASA DFRC contains the 
five modules described immediately below. 

The input module is used to read the user input data, which includes model type, 
analysis type, constraints, responses, optimizer setting, et cetera. 

The optimizer module is used to select any of the optimizers described in this 
paragraph, based on user input. Optimizers include a genetic algorithm (GA) with 
a continuous design variable (GC); a gradient-based algorithm, that is, DOT, with a 
continuous design variable (DC); a genetic algorithm with a discrete design variable (GD); 
GC + DC; and GC + DC + GD. Additional optimizers can be added in this module if they 
are needed in the future. 

The pre-processor module is used to create and update input files based on the 
design variable values before executing the analyzer module. 
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The analyzer module is used to execute analyzer and scripts automatically. Users 
can use a script file to execute a series of analyses in sequential order. 

The post-processor module is used to post-process the analyzer output file and 
to analyze the results automatically. Since the output file format from MSC Nastran is 
different than that from ZAERO® or CFL3D, several sub-modules have been developed. 
More sub-modules will be added as needed. 

Structural Analysis Discipline 

The structural analysis discipline has been developed into the MDAO framework; 
other disciplines are planned to be included in this framework in the future. The structural 
analysis discipline determines the structural sizing required for a vehicle to be of sufficient 
strength to withstand the load conditions imposed during operation. The structural design 
of an airframe in MDAO is determined by multidisciplinary criteria which include stress, 
deflection, buckling, and margins of safety for sizing optimization under given loading 
conditions; flutter and divergence requirements; and weight. This is a complex effort, which 
involves loads analysis, structural dynamics, aeroelasticity, and structural optimization. 
The main outputs from the structural analysis discipline are structural weight, mass 
properties, safety factors, divergence speeds, flutter speeds, flutter frequencies, natural 
frequencies, and mode shapes. 

Several sub-modules were developed in the post-processor module for the structural 
analysis discipline. These sub-modules automate data calculation, provide data output in 
useful formats, and provide information with which to compute the objective function and 
penalty function. 

The Weight Calculation Sub-module 

Weight calculation by this module is a straightforward process. MSC Nastran with a 
weight generator is used to compute the model weight, moment of inertias, and location 
of the center of gravity (CG). The weight will be saved as an objective function if defined 
so by the user. 

The Deflection, Stress, Strain, and Buckling Sub-module 

Deflection, stress, and strain are computed using the MSC Nastran static analysis 
solution sequence 101. For overall linear buckling, the MSC Nastran solution sequence 
105 is used to calculate the eigenvalue associated with the first buckling mode for the 
structure. The critical buckling load can also be calculated for the structure. 
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The Modal Analysis Sub-module 

Natural frequencies and mode shapes are computed using MSC Nastran solution 
sequence 103. These data are generated and saved for the ZAERO® flutter analysis. 

The Flutter Sub-module 

For higher accuracy, subsonic and supersonic flutter analyses are performed using 
the ZAERO® code. Since the ZAERO® code requires structure modal data, the modal 
analysis sub-module must execute first. 

The Divergence Sub-module 

Divergence analysis is performed using an in-house code based on eigenvalue 
problems. The theory is based on the Rayleigh-Ritz Method and Galerkin’s Method 
(ref. 11). Input data to this code is a generalized aerodynamic influence coefficient (AIC) 
matrix from the ZAERO® code, and generalized stiffness matrix data from the MSC 
Nastran code. The output is the divergence speed. 

APPLICATIONS 

The MDAO tool has been validated by several test cases and applications, including 
the following two applications and structural model tuning (ref. 10). 

Three-Bar Truss Stress/Strain Analysis Model 

The preliminary application of the MDAO tool was the optimal design of a three- 
bar truss subjected to an external load (ref. 8) as shown in figure 3. In this problem, the 
objective is to minimize the total weight of the structure. In figure 3, the design variables 
X 1 and X 2 correspond to the cross-sectional areas of member 1 (and 3) and member 2, 
respectively. The area of member 3 is “linked” to be the same as member 1 for symmetry. 
The constraints are tensile and compressive stress constraints in member 1 and member 
2 under loading P 1 and P 2 . The loading P 1 and P 2 are applied separately and the material 
specific weight is 0.1 lb/in 3 . The allowable stress of the tension and compression in the 
member is 20,000 psi and -15,000 psi, respectively. This problem can be solved by using 
a finite element model and an analytical model with a closed-form solution. Table 1 shows 
the comparison of results between the closed-form solution, MSC Nastran using solution 
200, the current MDAO tool with the DOT optimizer, and the current MDAO tool with the 
GA optimizer. All of the results agree with the analytical results. The optimization history 
of the objective function using the GA optimizer is shown in figure 4. 
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Ikhana with the Fire Pod 


The three-bar truss is a straightforward application. A more interesting application and 
one of interest to NASA DFRC is the analysis and optimization of aeroelastic surfaces. 

Overview 

The NASA DFRC is NASA’s premier flight research center and is often the location for 
flights of highly modified vehicles intended to explore highly unusual flight regimes. When 
a vehicle has been modified elsewhere, it is often DFRC’s task to analyze its aeroelastic 
characteristics for flight safety validation , and the same is true when the aircraft is designed 
or modified “in-house.” Aircraft must be free from aeroelastic instabilities such as flutter 
and divergence to ensure safe operation. 

Recently, NASA DFRC acquired a Predator B (General Atomics Aeronautical 
Systems, Inc., San Diego, California) unmanned aircraft system for civilian missions. The 
Ikhana, as the vehicle is named, is shown in figures 5 and 6 and will carry a “fire pod” that 
will transmit images of remote areas of the western United States down from the aircraft 
to a ground station. The fire pod is located under the wing, near the left wing root, and can 
alter the flutter characteristics of the baseline aircraft. The flutter flight envelope prediction 
of the Ikhana using the current MDAO design process is the second optimization problem 
of this study. 

The objective of the second optimization problem is to maximize the flutter and 
divergence speeds of the structure by varying the chordwise location of the fire pod 
+/-30 in. chordwise from the baseline location while constraining the allowable tensile and 
compressive stresses within an acceptable range. The chordwise location of the fire pod 
will be the design variable for this application. A negative chordwise location means to 
move the fire pod forward from the baseline location. 

Below are the steps of the MDAO flutter and divergence optimization for the Ikhana 
problem. Several structural analysis disciplines will execute in sequential order within the 
MDAO analyzer module for this optimization problem. 

Step 1: update the finite element model and the aerodynamic model based on the 
design variable in the pre-processor module. 

Step 2: run the modal analysis using MSC Nastran solution 103 to compute stiffness, 
mass, frequency and mode shape in the analyzer module. 

Step 3: run the ZAERO® trim analysis and MSC Nastran solution 101 to compute the 
stresses in the analyzer module. 

Step 4: run the ZAERO® to compute flutter speed in the analyzer module. 
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Step 5: run the in-house divergence code to compute divergence speed in the 
analyzer module. 

Step 6: run the post-processor to collect and compute flutter and divergence speeds 
and constraints. 

Step 7: return the objective function, constraints, and penalty function to the 
optimizer. 

Step 8: repeat from step 1 until the optimum result is reached. 

Challenge 

A flutter analysis determines the dynamic stability of an aeroelastic system. As with 
static aeroelastic analysis, flutter analysis presupposes a structural model, an aerodynamic 
model, and their interconnection by splines. Therefore, modification of the fire pod location 
affects both the structural finite element model and the unsteady aerodynamic model in 
this design optimization process. To achieve the true optimum result, new MSC Nastran 
and ZAERO® analyses must be executed for each optimization iteration with any design 
variable update. The challenge of the Ikhana example is in the size of the AIC matrices 
and the complexity of the problem. With current computing resources at NASA DFRC, 
and using the ZAERO® code to generate the AIC matrix, computation for a single case, 
one Mach number with 16 reduced frequencies, takes an average of 20 hr. Even using 
an existing AIC matrix, it still takes approximately 25 min. to complete the flutter analysis. 
Using the GAoptimizer requires at least hundreds and perhaps thousands of MSC Nastran 
and ZAERO® executions. This is not at all practical as a timely process. 

Approach 

As mentioned above, because of computing resource limitations and the excessive 
real time required for this large optimization problem, several approaches have been 
investigated to expedite the flutter and divergence calculations and to avoid computing a 
new AIC matrix for each design variable update. 

The first approach computes an AIC matrix that corresponds to the desired fire 
pod location using the cubic spline from some precomputed AIC matrixes based on a 
predefined fire pod location that has been previously investigated. A matrix named “AJJ” 
that contains the unsteady pressure coefficients (Cp) on each aerodynamic box for each 
Mach number and reduced frequency can be output from ZAERO®. Matrix AJJ is the 
fundamental matrix for AIC matrix calculation. Since matrix AJJ is very large, it is very 
time-consuming to cubic-spline each element in the matrix. In addition, ZAERO® does 
not allow users to modify or replace matrix AJJ. ZAERO® does allow the user to replace 
matrix QHH, which is the generalized aerodynamic forces matrix due to mode shapes. 
Since QHH is a generalized matrix, the QHH matrix is quite small compared to matrixAJJ. 
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Although matrix QHH can be calculated from matrix AJJ (ref. 3), computing matrix QHH 
from matrix AJJ requires much computing memory and power. The first approach was not 
chosen for the Ikhana problem due to constrained computing resources and excessive 
computational time. 

The second approach is to perform the approximation based on the matrix QHH. 
The decision was made to cubic-spline each element in the generalized aerodynamic 
forces matrix QHH (instead of AJJ) in the most efficient and effective way. The matrix 
QHH calculation was carried out explicitly for seven fire pod chordwise locations with one 
Mach number and 1 6 reduced frequencies. These generalized (modal) aerodynamic force 
coefficient matrices were then interpolated to any new fire pod chordwise location for the 
flutter analysis. The flutter analysis is performed in modal coordinates. The analysis does 
not provide accurate estimates of the AIC for the desired fire pod location. Using this 
second approach, the time required for computing the divergence speed and the flutter 
speed calculation using ZAERO® is still time-consuming, requiring 40 min. Because of the 
difficulty of using matrix AJJ or QHH, a third approach was developed and the decision 
was made to perform the approximation based on the third approach. 

The third approach, since the technical goal in this research is to demonstrate the 
aeroelastic optimization problem, is to determine the interpolated flutter and divergence 
speeds from some precalculated flutter and divergence speeds computed from ZAERO® 
and in-house codes based on some predefined fire pod locations without running MSC 
Nastran, ZAERO®, and in-house codes for each design variable update. 

Thirteen ZAERO® flutter and in-house divergence analyses with chordwise fire pod 
locations +/-30 in. with 5-in. increments from the baseline model were analyzed before 
optimization. The flutter and divergence speed are summarized and shown in figure 7 
and table 2 with case 7 being the baseline model. According to these results at thirteen 
chordwise pod locations, with 2% damping, in cases 9, 11, and 13, the first flutter speed 
decreased 57.7 knots equivalent airspeed (KEAS), 67.6 KEAS, and 173.1 KEAS while 
moving the fire pod 10 in. , 20 in., and 30 in. aft, respectively. In case 5, when the fire pod 
was moved 10 in. forward of the baseline, first the flutter mode became a hump mode 
then disappeared; divergence was the only instability. Based on figure 7 the divergence 
is critical for Ikhana with the fire pod between 25 and -30 in. This approach does provide 
fair estimates of the critical divergence speed for the desired fire pod location. This third 
approach enables performing the optimization using both the gradient-based optimizer 
and genetic algorithm. 


RESULTS AND DISCUSSION 

The Ikhana with fire pod optimization problem has been performed using both the 
gradient-based optimizer DOT and the genetic algorithm with the third approach. 
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Because of the complexity of this application and the presence of some local minima 
and maxima, the gradient-based optimizer works, although inefficiently, to reach the 
optimum solution for this problem. With different initial values, the optimizer can only 
reach different local maximum, not the global maximum. See table 3 for the results of the 
final design variable and the objective value using DOT. 

Multidisciplinary optimization of maximum flutter speed and divergence speed with 
no constraints using the genetic algorithm optimizer is exercised and the results have 
been presented in this report. The fire pod chordwise location is limited to -20 and 30 in. 
When the optimization was carried out with the third approach and the genetic algorithm 
optimizer with 20 populations, the optimum fire pod location was 14.4 in. aft of the baseline 
model and the divergence speed was 465.53 KEAS. Although there is no significant 
change in the value of divergence speed at each fire pod location, MDAO is able to 
find the optimum fire pod location in this exercise. Table 4 shows the comparison of the 
critical speed before and after optimization. The optimization result moved the V-g and 
V-u) curves from the baseline configuration in figures 8 and 9 to the optimized configuration 
in figures 10 and 11 . Only modes with frequency below 20 Hz are plotted. 

CONCLUSION 

A multidisciplinary design, analysis, and optimization (MDAO) tool using gradient- 
based optimizers and a genetic algorithm has been developed at the National Aeronautics 
and Space Administration Dryden Flight Research Center (DFRC) (Edwards, California) 
to automate the design and analysis process and leverage existing tools to enable true 
multidisciplinary optimization in the preliminary design stage of subsonic, transonic, 
supersonic, and hypersonic aircraft. The first stage of implementing and applying MDAO 
techniques to the fire-pod-equipped Ikhana has been successfully completed at DFRC; 
the MDAO shows much promise for improving the Ikhana-with-fire-pod optimization. 
Studies have also shown the potential of the MDAO tool for increasing the flutter margin 
of safety in flight-test envelope expansion. The MDAO process produces a much-desired 
design for flight envelope expansion. Applying the MDAO tool in modern and realistic 
aircraft design is efficient; once all of the tools for pre- and post-processing were in place, 
all design conditions were taken into account. 

The challenge of integrating higher-fidelity codes with the MDAO tool is in the time 
required for computation. High fidelity would become practical if high-fidelity analysis 
executed more quickly. All of the individual analysis components, as well as parts of the 
integrated procedure, have been numerically tested, with computing resource upgrades 
to improve computational time being carried out. Once the structural analysis discipline 
is finalized and integrated completely into the MDAO process, other disciplines such as 
aerodynamics and flight controls will be integrated. 
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TABLES 


Table 1 . Results of optimization on a three-bar truss. 


Bar area 

Closed-form 

solution 

MSC Nastran 

MDAO with 
DOT 

MDAO with 
GA 

Bar 1 

.788675 

.77142 

.78798 

.72093 

Bar 2 

.408249 

.45185 

.40999 

.66093 

Bar 3 

.788675 

.77142 

.78798 

.72093 

Total weight 

2.63896 

2.6338 

2.6388 

2.7003 

50 generation 

Number of iterations 

N/A 

5 

7 

and 

50 propulsion 


Table 2. Ikhana with fire pod divergence and flutter speed. 


Case 

number 

Fire pod 
location from 
baseline, in. 

Divergence 
speed, KEAS 

Flutter speed, 
KEAS 

Flutter 

frequency, Hz 

1 

-30 

466.5 

N/A 

N/A 

2 

-25 

465.2 

N/A 

N/A 

3 

-20 

464.0 

N/A 

N/A 

4 

-15 

463.2 

N/A 

N/A 

5 

-10 

463.3 

N/A 

N/A 

6 

-5 

463.7 

558.4 

14.6 

7 

0 

464.2 

554.2 

14.7 

8 

5 

464.5 

522.0 

14.6 

9 

10 

465.1 

496.5 

13.8 

10 

15 

466.5 

465.6 

13.3 

11 

20 

465.1 

486.6 

12.2 

12 

25 

464.7 

456.0 

10.8 

13 

30 

464.7 

381.1 

10.4 
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Table 3. Ikhana with fire pod optimization using design optimization tools. 


Design variable 

Chordwise location from baseline, in. 

Objective value 
Critical speed, KEAS 

Initial value 

Final value 

-20.0 

-19.80 

463.95 

0.0 

0.0001 

464.20 

10.0 

14.31 

465.53 

20.0 

14.27 

465.53 


Table 4. Summaries of critical speeds before and after optimization using 
genetic algorithm. 



Design variable 

Objective value 


Chordwise location 

Critical speed, 


from baseline, in. 

KEAS 

Before 

0.00 

464.20 

After 

14.37 

465.53 


11 


FIGURES 



Figure 1 . Central executive overview. 



Figure 2. The basic framework of the multidisciplinary design, analysis, and optimization 
tool central executive. 
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Figure 3. Three-bar truss load conditions. 
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Figure 4. History of design optimization of the three-bar truss: objective function weight 
and constraint stress using the genetic optimizer. 
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Figure 5. The Ikhana baseline aircraft. 



Figure 6. The Ikhana aircraft carrying the fire pod. 
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IKHANA with fire pod flutter and divergence speed 
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Figure 7. Ikhana with fire pod divergence and flutter speed. 


V-g plots 



Figure 8. The V-g curves before optimization. 
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V- co plots 



Figure 9. The V-oo curves before optimization. 


V-g plots 



Figure 10. The V-g curves after optimization. 
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Figure 11 . The V-co curves after optimization. 


17 


REFERENCES 


1 . MSC.NASTRAN 2005 Quick Reference Guide, MSC. Software Corporation, 2004. 

2. Matlab 7 Getting Started Guide, The MathWorks, Inc., 2008. 

3. ZAERO Version 8.0 User’s Manual, ZONA Technology, Inc., 2007. 

4. Krist, Sherri L., Robert T. Biedron, and Christopher L. Rumsey, CFL3D User’s Manual 
(Version 5.0), NASA-TM-1 998-208444, 1998. 

5. DOT Design Optimization Tools User’s Manual Version 5.0, Vanderplaats Research 
& Development, Inc., 2001. 

6. Vanderplaats, Garret N., Numerical Optimization Techniques for Engineering Design, 
with Applications, McGraw-Hill Book Company, 1984. 

7. Charbonneau, Paul, and Barry Knapp, A User’s guide to PIKAIA 1.0, 
NCAR/TN-418+IA, 1995. 

8. Moore, Gregory J., MSC/NASTRAN Design Sensitivity and Optimization User’s 
Guide, Version 67, The MacNeal-Schwendler Corporation, 1993. 

9. Goldberg, David E., Genetic Algorithms in Search, Optimization, and Machine 
Learning, Addison-Wesley Professional, 1989. 

10. Lung, Shun-fat, and Chan-gi Pak, Structural Model Tuning Capability in 
an Object-oriented Multidisciplinary Design, Analysis, and Optimization Tool, 
NASA/TM-2008-214640, 2008. 

11. Rao, S. S., The Finite Element Method in Engineering, Pergamon Press, 1982. 


18 


REPORT DOCUMENTATION PAGE 


Form Approved 
OMB No. 0704-0188 


The public reporting burden for this collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing 
data sources, gathering and maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or 
any other aspect of this collection of information, including suggestions for reducing this burden, to Department of Defense, Washington Headquarters Services, 
Directorate for Information Operations and Reports (0704-0188), 1215 Jefferson Davis Highway, Suite 1204, Arlington, VA 22202-4302. Respondents should be aware 
that notwithstanding any other provision of law, no person shall be subject to any penalty for failing to comply with a collection of information if it does not display a 
currently valid OMB control number. 

PLEASE DQ NQT RETURN YQUR FORM TQ THE ABOVE ADDRESS ¥ 


1. REPORT DATE (DD-MM-YYYY) 2. REPORT TYPE 

01-02-2009 Technical Memorandum 

3. DATES COVERED (From - To) 

4. TITLE AND SUBTITLE 

Multidisciplinary Design, Analysis, and Optimization Tool Development 
Using a Genetic Algorithm 

5a. CONTRACT NUMBER 

5b. GRANT NUMBER 

5c. PROGRAM ELEMENT NUMBER 

6. AUTHOR(S) 

Chan-gi Pak; Wesley Li 

5d. PROJECT NUMBER 

5e. TASK NUMBER 

5f. WORK UNIT NUMBER 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 
NASA Dryden Flight Research Center 
P.O.Box 273 

Edwards, California 93523-0273 

8. PERFORMING ORGANIZATION 
REPORT NUMBER 

H-2921 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 
National Aeronautics and Space Administration 
Washington, DC 20546-001 

10. SPONSORING/MONITOR'S ACRONYM(S) 

NASA 

11. SPONSORING/MONITORING 
REPORT NUMBER 

NASA/TM-2009-2 14645 


12. DISTRIBUTION/AVAILABILITY STATEMENT 

Unclassified -- Unlimited 

Subject Category 05 Availability: NASA CASI (301) 621-0390 Distribution: Standard 


13. SUPPLEMENTARY NOTES 

Pak, Li, NASA Dryden Flight Research Center. Also presented as ICAS08-951 at the 28th International Congress of hte Aeronautical Sciences (ICAS 
2008), Anchorage, Alaska, September 14-19, 2008. An electronic version can be found at http://dtrs.dfrc.nasa.gov or http://ntrs.nasa.gov/search.jsp 


14. ABSTRACT 

Multidisciplinary design, analysis, and optimization using a genetic algorithm is being developed at the National Aeronautics and 
Space Administration Dryden Flight Research Center (Edwards, California) to automate analysis and design process by leveraging 
existing tools to enable true multidisciplinary optimization in the preliminary design stage of subsonic, transonic, supersonic, and 
hypersonic aircraft. This is a promising technology, but faces many challenges in large-scale, real-world application. This report 
describes current approaches, recent results, and challenges for multidisciplinary design, analysis, and optimization as demonstrated 
by experience with the Ikhana fire pod design. 


15. SUBJECT TERMS 

Flutter, Genetic Algorithm, MDO, Multidisciplinary, Optimization 


16. SECURITY CLASSIFICATION OF: 

17. LIMITATION OF 
ABSTRACT 

18. NUMBER 
OF 

19b. NAME OF RESPONSIBLE PERSON 

a. REPORT 

b. ABSTRACT 

c. THIS PAGE 


PAGES 

STI Help Desk (email: help(a)sti. nasa.gov) 






19b. TELEPHONE NUMBER (Include area code) 

u 

U 

uu 

U 

23 

(301)621-0390 


Standard Form 298 (Rev. 8-98) 

Prescribed by ANSI Std. Z39-18 

































